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Abstract 

The  foundation  of  the  nonlinear  theory  of  asymmetric  anisotropic  sandwich  plates  with  a  first 
order  compressible  weak  orthotropic  core  under  a  Friedlander-Type  explosive  blast  is  presented. 
The  equations  of  motion  are  developed  by  means  of  Hamilton’s  Principle.  Within  the  theory,  the 
face  sheets  are  asymmetric  while  adopting  the  Love-Kirchoff  model.  In  addition,  the  core  layer 
is  assumed  to  be  compressible  (extensible)  in  the  transverse  direction  thereby  capturing  any 
wrinkling(local)  or  global  instabilities.  The  theory  is  then  simplified  and  applied  for  applicable 
cases  of  sandwich  plates  with  symmetric  orthotropic  and  cross-ply  facings  under  blast  loading. 
The  governing  solution  is  developed  using  the  Extended-Galerkin  method  resulting  in  two 
coupled  nonlinear  second-order  ordinary  differential  equations  which  are  then  solved  using  the 
Adaptive  (variable  time  increments)  4th-Order  Runge-Kutta  Method  for  a  system  of  differential 
equations.  Results  are  then  validated  for  the  case  of  a  unifonn  pressure  pulse  for  the 
incompressible  core  case.  It  is  shown  that  reasonable  agreement  exists. 

Key  Words:  Sandwich  Panel;  Dynamic  Response;  Blast;  Transient  Response;  Friedlander;  Shock 
Loading;  Structural  Response 

1.  Introduction 

During  combat  situations,  the  structure  of  army  military  vehicles  may  have  to  withstand 
explosive  blast  loading,  the  residual  effects  of  which  are  the  blast  pressure,  fragmentation,  and 
heat.  Due  to  their  outstanding  characteristics,  sandwich  structures  can  be  adopted  in  armored 
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plating  and  the  vehicle  hull  structure  of  armored  vehicles  and  possibly  other  key  parts  of  the 
vehicle.  Their  advantages  consist  of,  among  other  things,  a  high  bending  stiffness  and  strength  to 
weight  ratio  and  are  lightweight  in  structure.  These  structures  also  provide:  a)  excellent  thermal 
and  sound  insulation;  b)  increased  durability  under  a  thermo-mechanical  loading  environment; 
and  c)  tight  thennal  distortion  tolerances.  As  a  result,  these  structures  may  reveal  many 
advantages  when  implemented  into  military  armored  military  vehicles. 

The  purpose  of  this  paper  is  not  only  for  the  benefit  of  the  Anny,  in  regards  to  structural 
applications,  but  also  to  fill  in  some  of  the  gaps  existing  within  the  literature.  Several  papers 
currently  existing  present  results,  regarding  the  blast  of  sandwich  panels.  Among  these  papers 
several  factors  are  considered  while  others  are  neglected.  The  following  authors  [3,4,10,11] 
considered  blast  loading  of  sandwich  panels  while  neglecting  the  nonnal  or  transverse 
deformation  of  the  core.  Authors  such  as  [5,6,7,12]  have  considered  the  defonnation  of  the  core 
for  both  under  water  and  in-air  explosive-type  loading  with  the  in-air  case  limited  to  a  sonic 
boom  and  triangular  pulse-type  loading.  In  addition,  the  theory  was  limited  to  symmetric 
sandwich  panels.  To  advance  beyond  this,  the  authors  [14,15]  considered  the  effect  of  a  higher- 
order  compressible  core  (2nd-order)  with  an  exponential  decay.  Within  these  authors  works  a 
different  theoretical  approach  was  taken  as  compared  with  [5,6,7,12],  Numerical  or  experimental 
results  [2,9,13,17]  have  shown  that  the  core  experiences  significant  deformation  during  sudden 
impact.  Therefore,  the  consideration  of  the  deformation  of  the  core  within  the  governing  theory 
is  imperative.  Finally,  to  capture  a  more  representative  model  for  in-air  blast  loading,  The 
Friedlander-type  model  should  be  utilized.  To-date  the  most  comprehensive  analytical  theory 
presented  on  sandwich  panels  is  given  by  [5,6,7,12],  In  few  of  these  papers  [6,7,12]  a  nonlinear 
theory  for  symmetric  anisotropic  sandwich  plates  and  shells  incorporating  a  compressible  core 
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exposed  to  under  water  blast  and  in-air  blast  loading  has  been  presented  with  the  above 
mentioned  limitations. 

This  paper  expands  beyond  the  previous  authors  by  presenting  and  laying  the  foundation  for  a 
geometrically  nonlinear  theory  of  asymmetric  sandwich  plates  with  a  first  order  compressible 
core  exposed  to  an  in-air  Friedlander-type  explosive  loading.  The  results  are  then  presented  for  a 
simplification  of  the  governing  equations  for  the  case  of  symmetric  orthotropic  and  cross-ply 
facings. 

2.  Basic  Assumptions  and  Preliminaries 

Presented  in  Fig.  1  is  a  geometrical  representation  of  a  sandwich  structure.  The  structure  consists 
of  asymmetric  facings  with  a  thick  core  between  the  facings.  The  face  sheets  are  assumed  to  be 
asymmetric  with  respect  to  the  global  mid-surface  (mid-surface  of  the  core).  The  face  sheet 

thicknesses  for  the  top  and  bottom  face  sheets  is  denote  by  t'j  and  tbf ,  respectively,  while  tc 
denotes  the  thickness  of  the  core.  A  Cartesian  coordinate  system  is  placed  as  shown,  in  Fig.  1, 
with  the  transverse  direction,  x3  assumed  positive  in  the  downward  direction. 

The  tangential  stiffness  of  the  sandwich  plate  is  assumed  large  which  implies  small  tangential 
deformations.  In  contrast,  large  defonnations  can  occur  in  the  transverse  direction.  The  face 
sheets  are  assumed  incompressible  while  the  core  is  assumed  to  be  compressible  in  the  transverse 
direction.  Consistent  with  the  concept  of  small  tangential  and  large  transverse  displacements,  the 
tangential  and  rotatory  inertias  are  neglected.  Where  as,  only  the  transverse  inertias  are 
considered.  Finally,  in  addition,  it  is  assumed  that: 

1 .  The  face  sheets  fulfill  the  Love-Kirchoff  assumptions  and  are  thin  compared  with  the 
core. 

2.  The  bonding  between  the  face  sheets  and  the  core  is  assumed  to  be  perfect. 
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3.  The  kinematic  boundary  conditions  at  the  interfaces  between  the  core  and  the  facings  are 
satisfied. 

4.  The  core  is  assumed  to  be  a  weak  ortho  tropic  transversely  compressible  core  carrying 
only  the  transverse  strains  and  the  nonnal  strain. 

5.  The  shock  wave  pressure  is  uniformly  distributed  on  the  front  face  of  the  sandwich  plate. 
3.  Kinematic  Equations 

3.1  Displacement  Field 

Consistent  with  standard  plate  and  shell  theory  and  adopting  the  Love-Kirchoff  assumptions 
the  displacement  field  for  the  top  and  bottom  facings  is  given  as: 

Top  Face  Sheets 


t  a  .  d 
Va  =  Ua  +Ua- 


tc+tf 


U3,a  - 


„  (c+tf  d 
X3  U3,a 


V3  ~u3  +  U3 

Bottom  Face  Sheets 


(lb) 


(2a) 


bad 

v3  =  u3  -u3  (2b) 

In  the  above  equations,  the  Greek  indices  have  the  range  1,  2,  while  the  Latin  indices  have  the 
range  1,  2,  3  and  unless  otherwise  stated,  Einstein’s  summation  convention  over  the  repeated 
indices  is  assumed.  Also,  (•) ,  denotes  partial  differentiation  with  respect  to  the  coordinates  xt , 

while  superscripts  t  and  b  indicate  the  association  with  the  top  and  bottom  facings  respectively. 

A  second  order  power  series  is  assumed  for  the  core  tangential  displacements  while  a  first 
order  polynomial  is  assumed  for  the  core  transverse  displacement  [5,6,7,12].  With  this  in  hand, 
the  core  displacement  field  is  represented  in  compact  notation  as 
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Z,  Z,  1 

>;=Yju'z'’  Bl='Zw?*‘ 


(3a-c) 


Satisfying  the  interfacial  continuity  conditions  at  the  interfaces  between  the  facings  and  the  core, 
namely: 

Top  face  sheet/core  interface  (x3  =  —tc/2) 


C  t  C  t 

Va=Va’  V3=V3 


Bottom  face  sheet/core  interface  (x3  =tc/ 2) 

c  b  c  b 
»«=»«>  V3=v3 

Results  in  the  following  displacement  field  for  the  core. 

d  2x, 


»a=Ua- 


f  t 

h  ^ 

f  t 

h  \ 

1 

„ ,  a 

tf  +  tf 

4 

U3,a 

4 

V 

y 

V 

J 

U3  ,a 


3-Ua  + 


(  J  tb  \ 

f 

,t  .b  A 

tf  4 -tf 

^  „ ,  a  . 

i 

2/ 

X3U3,a  + 

2/, 

V  c  ) 

L 

c  y 

X3U3,a 


O' 


(4a, b) 


(5a,  b) 


(6a) 


v%(x,y,z,t)  =  u%  - 


2x 


3  ..d 


(6b) 


For  the  special  case  of  symmetric  facings  with  respect  to  the  global  mid-surface,  t‘f  =  tfj  =  t 1  .  In 


Equations  (la,b),  (2a, b),  and  (6a, b)  the  displacement  functions 


u i  =  ^ (ui  +  )»  Ui  =\ ^  ~U‘  ^ 


(7a, b) 


represent  the  average  and  the  half  difference  of  the  face  sheet  mid-surface  displacements  while, 
the  core  displacements,  O',  represent  warping  functions  of  the  core. 

3.2  Non-Linear  Strain-Displacement  Relationships 
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The  strain-displacement  relationships  given  by  the  Lagrangian  Strain-Displacement 
Relationships  used  in  conjunction  with  the  Von-Karman  assumptions  is  given  in  compact 
indicial  notation  as: 

27ij  =  vij  +  vjj  +  v3jv3j  (8a) 


In  expanded  form,  they  are  given  as 


y it  -  uu  +  2  (°3,i) 

(8b) 

1  ,  ,2 

722  _  V2,2  +  ^  (y3,2) 

(8c) 

1  ,  ,2 

733  _  v3,3  +  2  (°3,3) 

(8d) 

1,  ^  1 

723  -  2  +  °3’2  ^  +  2  U3>2°3’3 

(8e) 

1,  ,  1 

7l3  -  2  ^L3  +  °3,l)  +  2  °3,1°3,3 

(81) 

1,  ^  1 

7l2  “  2  (UL2  +  °2,l)  +  2  w3,ly3,2 

(8g) 

Substituting  the  displacement  relationships,  Eqs.  (la,b),  (2a, b)  and  (6a, b)  into  the  nonlinear 
strain-displacement  relationships,  for  the  individual  layers  of  the  structure  (top,  bottom,  and  core 
layers),  results  in: 


t  —a  ,  —d  . 

( 

tc+tf] 

a 

f 

tc+t'A 

lap  -  yap  +  yap  + 

x3  +■ 
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(9a) 


and 
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b  —a  —d  . 

yap  -  Yap  ~  yap  + 
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v  _ 

*c  +  fbf  ^ 

a 

Kap  ~ 

( 

V  - 

*c  +  ff 

x3 

V 

2 

9 

X3 

V 

2 

J 
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afi 


(9b) 


Where 


7  op  -  2  +  yafi )’  VaP  ~  2^°P  1 V )  (10a) 

and 

K°ap  =  ^(K'ap  +  ^apl  Ktp  =  (10b) 


In  the  above  expressions,  are  referred  to  as  the  average  and  half  difference  of  tangential  or 
membrane  strains  of  the  top  and  bottom  facings;  while,  K(^f  )  are  referred  to  as  the  average  and 

half  difference  of  the  bending  strains  of  the  top  and  bottom  facings.  The  expressions  for  the 
membrane  and  bending  strains  are  provided  in  Appendix  A. 

For  the  core,  the  strain-displacement  relationships  take  the  form 


7a 


—c,  c 
'  7/3  +ZK 


i  3 


(ID 


In  these  expressions,  y?3  and  xci3  are  the  membrane  and  bending  strains,  respectively.  These 
expressions  are  also  provided  in  Appendix  A. 

4.  Constitutive  Equations 

Both  the  top  and  bottom  face  sheets  are  considered  to  be  constructed  from  unidirectional  fiber 
reinforced  anisotropic  laminated  composites,  the  axes  of  orthotropy  not  necessarily  being 
coincident  with  the  geometrical  axes  (xt ,  x2  ) .  The  stress-strain  relationships  for  each  lamina  of 
the  facings  becomes 
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Tn  611  Qn  Q\6  7  ii 

x22>~  Q22  Q26  '  722 

T12  Sym  Q66  2y12 


Where 

Q  -  for  i,  j  =  (1,  2,  6)  are  the  Transformed  plane-stress  reduced  stiffness  measures.  See  Reddy 


The  stress-strain  relationships  for  the  orthotropic  core  with  the  geometrical  and  material  axes 
coincident  are  expressed  as 


^=Ecy  3C3,  r2c3  =  G2c3y2c3 . 


(13a-c) 


5.  Equations  of  Motion  and  Boundary  Conditions 

An  energy  approach  is  taken  to  determine  the  equations  of  motion  and  as  a  by  product,  the 
boundary  conditions.  Letting  U  represent  the  strain  energy,  W represent  the  work  done  by 
external  forces,  and  T represent  the  kinetic  energy,  Hamilton’s  Variational  Principle  can  be 


expressed  as 


fh 

(SU  -  SW  -  ST)dt  =  0 . 

J  tn 


It  should  be  mentioned  that  hamilton’s  principle  is  valid  for  both  elastic  and  elastic-plastic 
theoretical  models.  Assuming  a  weak  compressible  core,  the  variation  of  the  strain  energy  and 
work  can  be  written  as 


m-tc/ 2  r+‘c/2  eh 

,  t^apSyapdx  3+I  r-3dyf3dx  3+| 


f  ^ybapdx 3  dA 


Where  r;/-  are  the  tensorial  components  of  the  second  Piola-Kirchoff  stress  tensor,  while  A  is 


attributed  to  the  area  of  the  sandwich  plate.  With  respect  to  the  work  done  by  external  loads,  only 
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the  external  work  done  by  the  transverse  loading  due  to  the  explosive  pressure  pulse  as  well  as 
energy  loss  due  to  structural  damping  will  be  considered.  All  work  due  to  tangential  loadings 
will  be  neglected.  The  SW  takes  the  form 


SW  =  ^{^2{xl,x2,t)dv\  +  qb(xl,x2,t)Svb  —  2C( b\dv\  -2Cci)\dv\  ~2CbvbSvb^jdA  (16) 

where  q‘(xl,x2,t)  denotes  the  transverse  pressure  loading  from  a  spherical  air-blast  and  C  is  the 
structural  damping  coefficient  per  unit  area  of  the  plate.  Discarding  tangential  and  rotatory 
inertia  effects,  the  variation  of  the  kinetic  energy  can  be  written  as 


r‘\  ft\  r  f  Md2  ,  ,  ,  r'c/2  r^/2+4  ,  .  , 

I  STdt  =  I  I  -  I  P  fvidvidx-,  +  1  pcv?dvZdx-,  +1  phfv^5v\dx 

4  JJa  (  J-tc/2-/f  f  3  J  -tc/2  ~  4/2  y 


dAdt 


(17) 


Where  pc  and  prf,  pbf ,  are  the  mass  densities  of  the  core  and  the  top  and  bottom  face  sheets, 
respectively,  and  v  denotes  the  transverse  acceleration. 

The  equations  of  motion,  in  terms  of  the  global  stress  resultants  and  stress  couples  are  obtained 
by  substituting  Eqs.  (15)  -  (17)  and  the  strain-displacement  relationships,  Eqs.  (9)  -  (1 1),  while 
utilizing  the  definitions  of  the  stress  resultants  and  stress  couples  into  Hamilton’s  Equation  (14), 
while  collecting  the  coefficients  of  each  of  the  virtual  displacements,  integrating  by  parts  where 
ever  feasible,  and  invoking  the  independent  and  arbitrary  character  of  the  virtual  displacements 
results  in  the  following  nonlinear  equations  of  motion  along  with  the  corresponding  boundary 


conditions. 

K  : 

=  0 

(18a) 

duda  : 

Nc 

A [d  ,  Y«3  _  A 

(18b) 

Ki  =  o 

(18c) 
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$u3  '■  u3,apNap  +Map  ap  +U3  a/3Naj3  + 
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U-,  - 


c‘  +ch 
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,d  ,  ^3+^3 


|m3  + 
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•  d 

«3  - 


-CM 
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In  the  above  equations,  the  global  stress  resultants  and  stress  couples  are  defined  as 


(18e) 


fc ,  )  = =  \  fc  +  N%  )  (m^  +  )} 

(< ,  M  J  )  =  i  -  <  )  (mJ,  -  )} 


(19a) 

(19b) 


Where  the  local  stress  resultants  and  stress  couples  are  given  as: 


(  )  p— fc/2 

J-trU-tf 


b+t'P 


J  J 


> dx , 


(20a) 


tr 


ted 


Xa/3  3 


X-> 


‘c+‘d 1 


(20b) 


f  'c/24(1,x3)^3,  (/  =  1,  2,  3) 

J-C/2 


(20c) 
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The  local  stress  resultants  and  stress  couples  can  be  expressed  in  terms  of  the  displacements  by 
substituting  the  corresponding  constitutive  relationships,  Eqs.  (12),  (13a-c)  along  with  the  strain- 
displacement  relationships,  Eqs.  (9)  -  (1 1)  into  Eqs.  (14)  -  (16).  This  results  in 
Top  face: 


Bottom  face: 
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For  the  special  case  of  symmetric  facings  with  respect  to  the  local  mid-surfaces,  =  0  .  The 


stress  resultants  in  terms  of  the  strain  deformation  components  of  the  core  take  the  fonn 
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(24a,  b) 


The  core  stiffness’s  are  given  by 

(A‘0,D‘0)=  rc/2Qfa(!,xi)dx3, 

J-tn/2 


i  =  ( 3,  4,  5) 


(25) 


Where  no  summation,  with  respect  to  the  repeated  indices,  is  assumed.  The  core  transverse  and 
nonnal  moduli  are  given  as: 

Qx=Ec,  Qc44=gcvz,  Qc55=gcxz 


(26) 


The  corresponding  boundary  conditions  are  given  as: 
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=  0 

<n=  0  or  Mann=  0 

u(n=0  or  A/'  =0 


(27g) 

(27h) 


Where  «  and  t  are  the  nonnal  and  tangential  directions  to  the  boundary.  When  n  —  1 ,  t  =  2,  and 
when  n  —  2  ,  t  =  1  .  The  quantities  with  an  overcaret  represent  the  prescribed  quantities  on  the 
boundary. 

For  the  case  of  simply  supported  boundary  conditions,  the  boundary  conditions  become: 
Along  the  edges  xn  =  (0,  Ln ) 


Kn  =  <  K  =  K  =  Ml  =  Mdnn  =  u“  =  ui  =  0  (28a-h) 

6.  Spherical  in-Air  blast  Loading  of  the  Friedlander-Type 

With  the  ever  increasing  demands  for  increased  safety  for  the  warfighter  in  the  field  to  operate 
structurally  sound  vehicles  in  the  event  of  an  IED  or  some  other  type  of  explosive,  it  is 
imperative  that  an  understanding  of  the  structural  response  of  various  components  within 
military  combat  vehicles  under  an  explosive  blast  be  understood  so  that  measures  can  be  taken 
from  a  design  standpoint  to  ensure  the  durability  and  survivability  of  these  components.  To  begin 
to  achieve  this  understanding,  the  type  of  explosive  loading  considered  here  is  a  free  in-air 
spherical  air  burst.  Such  an  explosion  creates  a  spherical  shock  wave  which  travels  radially 
outward  in  all  directions  with  diminishing  velocity.  The  fonn  of  the  incident  blast  wave  from  a 
spherical  charge  is  shown  in  Fig.  2.  Where  qso  is  the  peak  overpressure  above  ambient  pressure, 
q0  is  the  ambient  pressure,  ta  is  the  time  of  arrival,  tp  is  the  positive  phase  duration  of  the  blast 
wave,  and  t  is  the  time.  The  waveform  shown  in  Fig.  2.  is  given  by  an  expression  known  as  the 
Friedlander  equation  and  is  given  as 
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q,  (0  =  (Vso  -  4o )[!  -  (t  -  ta )  /  tp  ]  exp [-«(f  -  ta  ) /  tp  ]  (29a) 

Where, 

g50  =  1772/Z3-114/Z2+108/Z  (29b) 

1/3 

In  Eq.  (29b),  Z  is  known  as  the  scaled  distance  given  by  z  —  R/W  with  R  being  the  standoff 
distance  in  meters  and  W  being  the  equivalent  charge  weight  of  TNT  in  terms  of  kilograms. 
Also,  a  is  known  as  the  decay  parameter  which  is  detennined  by  adjustment  to  a  pressure  curve 
from  a  blast  test. 

For  the  conditions  of  standard  temperature  and  pressure  (STP)  at  sea  level,  the  time  of  arrival  ta 
and  the  positive  phase  duration  tp  can  be  determined  from  [8] 

t/tx  =R/Rl  =(W/Wl)113  (30) 

Where  t\  represents  either  the  arrival  time  or  positive  phase  duration  for  a  reference  explosion  of 
charge  weight  W\ ,  and  t  represents  either  the  arrival  time  or  positive  phase  duration  for  any 
explosion  of  charge  weight  W.  The  detennination  of  the  standoff  distance  for  any  charge  weight 
W  follows  a  similar  reasoning.  The  application  of  these  relationships  is  known  as  cube  root 
scaling.  It  should  be  understood  that  in  applying  these  relationships  that  the  standoff  distances 
are  themselves  scaled  according  to  the  cube  root  law. 

7,  Solution  Methodology 

Up  to  this  point,  the  overall  governing  system  of  equations  developed,  applies  to  asymmetric 
sandwich  plates  both  locally  and  globally  with  laminated  composite  facings  and  a  weak 
transversely  orthotropic  first-order  compressible  core  under  transverse  loading.  For  the  purpose 
or  scope  of  this  paper,  the  equations  will  be  simplified  for  the  case  of  symmetric  sandwich  plates 
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with  orthotropic  and  cross-ply  facings  with  a  transversely  orthotropic  compressible  core  under  a 
Friedlander-type  of  explosive  transverse  loading. 

6.1  Special  Case:  Symmetric  orthotropic  single  layer  facings 
For  this  special  case,  t'j-  =  thf  =  t  f  .  In  fulfillment  of  the  geometric  boundary  conditions,  a 

suitable  representation  for  if ,  and  if  is  given  by: 


«3  =  wmn  (0  sin(4x! )  sin  (jinx2 )  (31a) 

w  3  =  WL  (0  sin(4xi )  sin(/<„x2  )  (31b) 

where  Xm  =  nm/ L] ,  pn  =  nn/  L2  .  m  and  n  are  the  number  of  sine  half-waves  in  the  corresponding 

directions  whereas  wamn  (t)  and  \ifin  (t)  denote  the  modal  amplitudes  as  a  function  of  time  of  the 
transverse  displacement  functions.  The  transverse  loading  is  represented  by 

qt(X],x2,t)  =  q  mn  (0  sin(/LwXj )  sin(yunx2  ),  (32a) 

which  implies  through  integration  of  both  sides  over  the  plate  area  that 

4  rL 2  rL\ 

q,»n  (0  =  ~r—  qt  (^1  ,X2  ,t)  sin(Amx, )  sin(/<„x2  )dxldx2  .  (32b) 

Jo  Jo 


Letting, 

qt  (x,  ,x2,t)  =  q,  (t)  =  ( qso  -  q0 )[1  -  (t  -  ta ) / tp ] e>q> [~a(t  -ta)/tp\ 
and  integrating  gives 


qm,tt) 


1 6q,(t) 

mnn 2 


(32c) 


(32d) 


The  chosen  solution  methodology  of  the  governing  system  of  equations  is  the  Extended 
Galerkin  method.  With  this  in  mind,  the  first  two  equations  of  motion,  Eqs.  (18a),  can  be 
satisfied  by  assuming  a  stress  potential  approach  with  the  following  representations. 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


N afi  ^'aoA' fjp(p  (jjp 


(33) 


where  cap  denotes  the  2-D  pennutation  symbol.  This  introduces  another  variable  in  the 

governing  system  of  equations  that  requires  determination.  The  detennination  comes  from  a 
compatibility  condition  or  equation.  The  compatibility  equation  can  be  obtained  from  the 
membrane  or  tangential  strains  from  the  strain-displacement  relationships,  Eqs.  (9,10)  and 
Appendix  A  by  eliminating  the  in-plane  displacements.  This  results  in: 

7l  1,22  —2y1212  +722,11  —  (w3,12)  +  (M3,12)  —  w3,l  lw3,22  —  w3,l  lw3,22  (34) 


Keeping  in  mind  that  for  symmetric  ortho  tropic  and  or  cross-ply  facings  \B]j'h)~\  =  0  and 


A  j/() ,  D/i6  =  0 ,  the  following  global  constitutive  relationships  can  be  determined  by  utilizing 


Eqs.  (19a,b)  and  (21a,b)  which  are  given  as 
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(35a, b) 


where  M  =  =  Ahu  .  By  performing  a  matrix  inversion  in  Eq.  (35a),  utilizing  Eqs.  (33),  and 


substituting  into  Eq.  (34)  provides  the  following  nonlinear  partial  differential  equation  in  terms 
of  the  Airy’s  stress  potential. 

^  ^  ^  ^  ex  2  d  2  a  a  d  d 

A\\(P, 2222  +(2^12  +^66  Ml  122  +^22^,1111  =  (W3, 12  )  +(W3,12)  —  W3,l  1W3,22  —  W3,l  1W3,22  (36) 


where  [A*  ]  =  [A?  ]  1 .  From  Eq.  (36)  cp  can  be  determined  by  assuming  the  following  general 


functional  form  then  substitituting  into  Eq.  (36),  and  comparing  coefficients. 
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cp(xx,x2,t)  =  (wamn)2cx  cos(2 kmxx)  +  (wamn)2c2  cos(2/<;)x2 )  +  (wdw)2 C3  cos(22mx,)  + 
(^L)2C4COS(2M„X2) 

where  the  expressions  for,  C, ,  C2  ,C3 ,  and  C4  are  provided  in  Appendix  B.  With  the  Airy’s  Stress 

potential  function  in  hand  along  with  Eqs.  (34a, b),  (33),  the  first  two  relationships  of  Eq.  (35a), 
and  the  expressions  for  the  membrane  strain  components  from  Appendix  A,  yields  a  system  of 

two  coupled  inhomogeneous  partial  differential  equations  in  terms  of  u“  and  u"  .  By  assuming 

the  following  two  general  solutions  for  u“  and  u2  ,  shown  below  in  Eqs.  (38a-b)  and 
substituting  these  expressions  into  the  coupled  partial  differential  equations,  and  comparing 
coefficients  gives. 


=  («„)2  +(wi)2)o,x1  +(«„)2  +<C„)2)d2  sin(2A„,x,)- 

(<  <„f  +  <  >*L  f  A  sm<  )««< 2/ vv ) 


(38a) 


i,J(x,,x2,0  =  (<<J2+<wi)2)£l^+((<„)2+<»i)2Ksin<2V2)- 

(«,» )2  +  <>*».) 2  cos(21„,x, )  sin(2  fi„x2 ) 


(38b) 


Where  the  constants  Dx  -  D3  and  A,  -  A3  are  provided  in  Appendix  B. 

By  expressing  the  third  and  fourth  equations  of  motion,  Eqs.  (18b)  in  terms  of  displacements 
utilizing  Eqs.  (24)  and  (35b)  with  the  corresponding  strain  deformation  components  from 
Appendix  A  two  coupled  partial  differential  equations  result.  These  coupled  partial  differential 

equations  can  be  solved  simultaneously  by  assuming  ux  ,  and  iti  in  the  following  form. 

U\  (A ,x2,t)  =  wamnwdmnAx  sin(2zmx, )  +  wamnwdmnA2  sin(22wx1 ) cos(2/i„.r2 )  + 
wamnA3  cos(2mx1)sin(yu„x2) 
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ui (x,  ,x2,t)  =  wamnwdmnBx  sin(2/<;)x2 )  +  wamnwdmnB2  cos(22mx1 ) sin(2/<„x2 )  + 

KnB 3  sin(2mx1)cos(/4„x2) 

Then  substituting  these  displacement  representations  into  the  coupled  partial  differential 
equations  and  comparing  coefficients.  These  coefficients  ^  -  A3  and  Bt  -  B3  are  provided  in 
Appendix  B. 

Utilizing  the  core  deformation  components  k‘b  from  Appendix  A  in  conjunction  with  the 

constitutive  equations  (24b)  the  fifth  and  sixth  equations  of  motion  (18c)  can  be  expressed  in 
tenns  of  displacements.  The  result  is  two  simultaneous  algebraic  equations  allowing  for  the 

detennination  of  .  The  result  is  given  as 


tc  d  1  d 

=  — M3  ,a~~u 


3  u3,a 


(40) 


This  relationship  reduces  an  9  parameter  system  ( uf ,  uf , , and  (p)  to  a  7  parameter  system  ( 

uf ,  uf ,  and  <p  )  by  eliminating  the  variable  (h  j  from  the  governing  system.  At  this  point  the  first 

six  equations  of  motion  are  satisfied  and  all  displacement  functions  are  known.  These 
displacement  functions  Eqs.  (31a,b),  (37),  (38a, b),  (39a, b),  and  (40)  satisfy  all  of  the  transverse 
geometric  boundary  conditions  and  the  boundary  conditions  with  respect  to  the  stress  couples, 
Mfm  and  Mdn  .  The  remaining  boundary  conditions  with  respect  to  the  tangential  stress 

resultants,  Nfm ,  Nf’t ,  Ndn ,  Ndt  will  be  satisfied  in  an  integral  average  sense  through  application 
of  the  Extended-Galerkin  Method.  The  only  unknowns  at  this  point  are  the  modal  amplitudes, 
wamn  (t)  and  wfnn  (t)  which  are  determined  by  means  of  the  extended  Galerkin  Method,  (Hause  et 
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al  [2,3]).  By  expressing  the  last  two  unfulfilled  equations  of  motion,  Eqs.  (18d,e),  and  the 
unfulfilled  boundary  conditions  in  terms  of  displacements  through  the  use  of  Eqs.  (33),  (35a, b), 
(37),  (38a, b),  and  (39a, b)  and  retaining  these  expressions  in  Hamilton’s  Energy  Functional 
carrying  out  the  indicated  integrations  and  collecting  the  coefficients  of  the  modal  amplitudes 

w“nn  and  wdmn  keeping  in  mind  that  the  variations  of  Swamn  and  Swdin  are  arbitrary  and 


independent  from  each  other  and  that  the  corresponding  coefficients  must  vanish  results  in  two 
nonlinear  coupled  second  order  ordinary  differential  equations  in  terms  of  the  modal  amplitudes 
which  is  solved  using  the  Adaptive  4th-order  Runge-Kutta  method  for  a  system  of  differential 
equations.  This  sytem  of  differential  equations  is  given  as: 


m,wa  +Cw“  +CcLwa  +  Cl\wa  wd  +  C?7wa  ( wd  )2  +  CUwa  Y  = 

1  mn  mn  10  mn  11  mn  mn  12  mn\  mn;  30  \  mn; 


,,a  \  3  9  r, 


(41a) 


WlWmn  +  C™mn  +  CqIWL  +  C02  (AW  +  Cq3  (%„,  )3  +  C$ 0  «„  f  +  C(x  «„  f  Wmn  =  (41b) 


The  coefficients  C10  -  C12 ,  C30 ,  C01  -  C03 ,  C20 ,  C21  are  expressions  which  depend  on  the  material 
and  geometrical  properties  of  the  structure  which  are  provided  in  Appendix  C.  C  is  the  damping 
coefficient,  while  qmn  is  the  amplitude  of  the  transverse  loading. 

8,  Results  and  Discussion 

Results  are  presented  next  for  the  dynamic  response  of  simply  supported  sandwich  plates  along 
all  four  edges  and  freely  movable  under  various  fixed  or  varied  geometrical  and  material 
properties  of  the  structure  exposed  to  blast  loading  to  develop  a  deeper  understanding  of  what 
effect  these  parameters  play  in  the  behavior  of  the  structure.  The  results  consist  of  two  parts.  The 
first  part  is  concerned  with  comparisons/validations  and  the  second  part  is  concerned  with  the 
current  or  present  results  produced  from  the  governing  theory. 

8.1  Validation 
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To  validate  the  results.  The  dynamic  response  of  a  sandwich  plate  with  an  incompressible  weak 
core  impacted  by  a  uniform  pressure  pulse  was  chosen  from  R.S.  Alwar  et  al  [1].  The  sandwich 
plate  is  assumed  to  be  simply  supported  and  freely  movable  along  all  four  edges.  The  uniform 

pressure  pulse  is  given  by  qmn(t)  =  \6q{) / tz2 mn  where  q0  is  the  unifonn  pressure. 

The  geometrical  and  material  properties  used  for  the  validation  are  given  for  the  face-sheets 
as,  Lx  =  L2  =  1.83  m,  tj-  =  0.000406  m , E(  =  68.96  (GPa),  E2  =  68.96  (GPa), 

G{2  =  25.92  (GPa),  v{2  =  0.33 ,  and  p  f  =  2768.38  (Kg/m3) .  The  corresponding  properties  for 
the  core  are  provided  as  tc  =  0.00635  m ,  GjC3  =  22.5914  (MPa),  G23  =  22.5914  (MPa), 
pc  =121.85  (Kg/m3)  .  Other  parameters  such  as  the  fiber  orientation  and  layup  of  the  facings, 
the  nondimensional  damping  parameter,  and  the  uniform  pressure  are  given,  respectively,  as 
[0/ Core/ 0] ,  A  =  0 ,  q0  =  61.6  Pa  and  123.2  Pa  .  From  Fig.  3,  it  can  be  seen  that  close  agreement 

exists.  It  should  also  be  noted  that  for  the  incompressible  core  case  wfnn  =  0  which  implies  that 

w“nn  (0  =  wmn(t)  which  is  constant  throughout  the  thickness  of  the  structure. 

8.2  Present  Results 

For  the  present  results,  unless  stated  otherwise,  the  following  geometrical  and  material 
properties  are  given,  respectively,  for  the  face-sheets  as  t  f  =  5mm ,  Lx  —  L2  -  900mm , 

E{  =  207  (GPa) ,  E{  =5.17  (GPa),  g{2  =2.55(GPa),  v{2  =0.25,  and  pf  =  l,588.22(K^m3) . 

The  corresponding  properties  for  the  core  are  provided  as  tc  =  50mm ,  E3  =3  (GPa) , 

G/3  =  0.1027  (GPa),  G23  =  0.0621  (GPa) ,  pc  =16  (KgTn3)  .  Other  parameters  such  as  the  fiber 
orientation  and  layup  of  the  facings,  the  nondimensional  damping  parameter,  the  peak 
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overpressure  qs0 ,  the  ambient  pressure  q0 ,  the  positive  phase  duration  tp  ,  and  the  rate  of  decay 
a  are  given,  respectively,  as  [O/Core/O],  A  =  0.05,  qs0  =1.766(MPa),  q0  =101(KPa), 
tp  =  0.001372(sec) ,  and  a  =  0.5  .  It  is  also  assumed  that  the  time  of  arrival,  ta  =  0  or  that  the 
time  of  the  blast  begins  at  t  —  0 . 

The  present  results  begin  with  Fig.  4  and  5.  depicting  the  effect  of  the  transverse  modulus  of 
the  core  on  the  global  and  wrinkling  response  of  the  sandwich  panel.  It  should  be  clarified  that 
the  terms  global  and  wrinkling  response  will  be  refered  to  throughout  the  analysis  of  the  results. 
The  global  response  is  a  measure  given  by  the  average  of  the  top  and  bottom  face  sheet 

deflection-time  histories  and  is  measured  by  w“m  (/) .  While,  the  wrinkling  (local)  response  gives 
a  measure  of  half  of  the  difference  between  the  top  and  bottom  face  sheets  with  regards  to  the 
deflection-time  histories,  w^m  ( t )  .  It  should  be  noted  that  in  the  case  of  the  incompressible  core, 

wfnn  (0  =  0.  Also,  in  sandwich  structural  applications  the  goal  from  a  design  standpoint  is  to 
contain  or  limit  the  amount  of  wrinkling  that  occurs  within  the  structure.  With  this  in  mind,  It  is 
revealed,  in  Figs.  4  and  5,  that  for  larger  values  of  the  young’s  modulus  of  the  core  the  effect  on 
the  global  response  is  minimal.  Whereas,  for  the  case  of  the  wrinkling  response  in  Fig.  5,  the 
effects  upon  the  deflection-time  histories  appears  to  be  more  beneficial  as  a  result  of  the 
decreased  amplitudes  of  deflection. 

In  Fig.  6,  it  can  be  seen  that  for  larger  rates  of  decay  that  the  global  deflection-time  history 
response  remains  for  longer  periods  of  time  within  the  negative  phase  or  suction  phase  of  the 
blast  pulse.  In  Fig  7,  which  is  the  counterpart  of  Fig.  6  for  the  case  of  the  wrinkling  response  it  is 
shown  that  a  similar  trend  is  observed  as  in  Fig.  6.  In  addition,  it  appears  that  the  frequencies  are 
larger  for  both  cases  as  the  rate  of  decay  is  increased. 
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Fig.  8  shows  that  global  deflection-time  response  is  enhanced  with  the  thicker  core  from  the 
standpoint  of  decreased  amplitudes  of  deflection.  Also  the  frequency  of  the  deflection-time 
response  begins  to  grow  out  of  phase  with  the  thinner  core  revealing  a  larger  frequency  of 
oscillation.  The  same  behavior  can  be  seen  in  Fig.  9  which  is  the  counterpart  of  Fig.  8  concerned 
with  the  wrinkling  response. 

Fig.  10  presents  the  effect  of  the  stacking  sequences  and  the  increase  in  thickness  of  the 
facings.  It  can  be  seen  that  the  cross-ply  layup  coupled  with  an  increase  in  the  face  thickness 
greatly  enhances  the  structural  response.  It  is  also  observed  that  the  cross-ply  layup  with 
increased  thickness  is  a  little  out  of  phase  with  the  orthotropic  single  layered  facings.  The  same 
analysis  also  holds  true  for  the  wrinkling  response.  It  appears  that  for  the  cross-ply  layup  with 
increased  face  thickness,  that  the  wrinkling  response  is  only  marginally  or  negligibly  affected. 

In  Figs.  12  and  13,  it  can  be  seen  that  the  global  response  is  more  sensitive  to  the  core 
transverse  modulus  ratio  than  the  wrinkling  response  in  Fig  13. 

9.  Concluding  Remarks 

The  governing  theory  of  asymmetric  sandwich  plates  with  a  first-order  compressible  core 
impacted  by  a  Friedlander-type  of  blast  has  been  presented  and  simplified  for  the  case  of 
symmetric  cross-ply  and  single-layered  orthotropic  facings.  In  all  cases,  it  was  mentioned  that  all 
four  edges  are  simply  supported  and  freely  movable.  Results  were  then  presented  for  this 
simplified  case  and  validated  against  results  found  in  the  literature  from  R.  S.  Alwar  et  al.  [1],  It 
was  found  that  for  the  incompressible  core  case  that  there  was  close  agreement  among  the 
results.  In  regards  to  the  compressible  core  case,  no  appropriate  results  have  been  found  in  the 
literature  for  the  theory  presented  in  this  paper  for  the  simply  supported  case  with  all  edges  freely 
movable.  The  effect  of  a  number  of  important  geometrical  and  material  parameters  were 
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analyzed  with  conclusions  drawn.  Some  of  the  important  conclusions  were  that  wrinkling 
response  seems  to  be  diminished  as  the  young’s  modulus  of  the  core  is  increased.  The  same  is 
the  case  for  larger  rates  of  decay.  Also,  for  thicker  cores,  both  the  global  and  wrinkling  responses 
are  less  severe.  It  was  also  revealed  that  the  compressibility  of  the  core  has  only  a  marginal  effect 
upon  the  global  response  of  the  sandwich  plate.  Finally,  the  cross-ply  type  layup  when  compared 
with  single-layered  facings  seemed  to  have  a  large  effect  on  the  global  response  and  less  effect 
on  the  wrinkling  response. 

One  should  keep  in  mind  that  both  the  stress  and  strain  profiles  should  be  detennined  to 
detennine  possible  failure  of  the  structure. 
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Appendix  A.  Strain-Displacement  Components 

Strain-Displacement  components  for  the  face-sheets 


7n  =  uii  +  2^3,1)  +  ^(.u3,^ 

722  =u2,2  +^(u3,2)  +^(u3,2) 

7i2  =  2Wi-2  +  2W2,1  +  2M3,iW3,2  +  2W3,iW3,2 
7ll  =  «u  +  w3,lw3,l 
722  =  u2,2  +  m3,2w3,2 


—  u  u.  u.  uu.  u  u 

7 12  -  2M1,2  +^'M2,1  +  2M3aW3,2  +  2M3,2m3,1 


The  bending  strains  for  the  face  sheets 


*11 

~  M3,ll 

a 

a 

k22 

-  -m3,22 

a 

a 

Kn 

-  ~U2,\2 

d 

d 

KU 

~  M3,ll 

d 

d 

K22 

-  M3,22 

d 

K\2 

-  ~U3,12 

Strain-Displacement  components  for  the  Core 


7 13  = 


1  </ 
- Mi 

C 


2tc  + 
4  C 


f3.1 


'V  h 

V _ V 

4C 


d  1  ..d-.a 
n3 1  —  —  Hi,  u 31 


723  -  ' 


- M  i  + 


4C 


l3,2 


V _ V 

4  C 


a 


2  1/ 


3  w3,2 


733 


2  d  2  d  2 

—  m3  +— (m3  ) 

C  t: 


*13  =- 


““^3,1  2"  ^3^3,1  +~<1,1 

C  t]  t] 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


c  1  d 

K23  = - 11 


d  d 


„  „  4  c 

3,2  +  2  l<3-2  +  t  ®2 

t  t 

lc  Lc 


K33  ~  ^ 


Appendix  B.  Constants 


Ci  -  C3  - 


(ftA[Ai4-(4)2] 
32  4j24 


(4)2[44-(4n 


A2- 


Ct  —  C/i  — 


32(yM„)  y4: 


2  A 
22 


4  = 


V 


c  7 


44+4 


Bn 

2 


Bn  ^22  ^mMl  + 


yz 


2?,= 


2  /  ^GL 

4/C  2  +  t  ' 
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Appendix  C.  Elastic  Coefficients 
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Figure  Captions 

Fig.  1  A  depiction  of  an  asymmetric  sandwich  plate  under  an  in-air  Friedlander  explosion. 

Fig.  2  Incident  profile  of  the  pressure -time  relationship  for  a  Friedlander-type  blast. 

Fig.  3  The  nondimensional  global  deflection-time  response  for  the  case  of  an  incompressible 
core  impacted  by  a  uniform  pressure  pulse 

Fig.  4  The  effect  of  the  transverse  modulus  of  the  core  on  the  global  response  of  a  sandwich 
plate  with  orthotropic  facings. 

Fig.  5  The  counterpart  of  Fig.  4  for  the  wrinkling  response  of  a  sandwich  plate. 

Fig.  6  The  effect  of  the  rate-of-decay  parameter  on  the  global  response  of  a  sandwich  plate  with 
ortho  tropic  facings. 

Fig.  7  The  counterpart  of  Fig.  6  for  the  wrinkling  response. 

Fig.  8  The  effect  of  the  core  thickness  on  the  global  deflection-time  history  of  a  sandwich  plate 
with  orthotropic  facings. 

Fig.  9  The  counterpart  of  Fig.  8  for  the  wrinkling  response. 

Fig.  10  The  effect  of  the  stacking  sequence  of  the  facings  on  the  global  response  of  a  sandwich 
plate. 

Fig.  1 1  The  counterpart  of  Fig.  12  for  the  wrinkling  response. 

Fig.  12  The  effect  of  the  core  shear  modulus  ratio  on  the  deflection-time  history  of  cross-ply 
laminated  sandwich  plate. 

Fig.  13  The  counterpart  of  Fig.  12  for  the  wrinkling  response. 
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Fig.  2 
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